AD  71X436 


Final  Report  on 

Quenching  of  Solid  Propellant  Combustion 
Grant  AF-AFOSR-0897-67 


March  1970 

Department  of  Chemical  Engineering 

This  research  under  Grant  AF-AFOSR-0897-67,  Project  task  No.  9711-01 
and  for  the  period  March  1967  through  February  1970  was  sponsored  by  the 
Air  Force  Office  of  Scientific  Research,  United  States  Air  Force. 

Technical  supervisor  for  this  program  is  Lt.  Col.  R.  W.  Haffner, 
project  scientist.  Propulsion  Division,  Directorate  of  Engineering 
Sciences,  Air  Force  Office  of  Scientific  Research. 

Reproduction,  translation,  publication,  use  and  disposal  in  whole  or 
in  part  or  for  the  United  States  Government  is  permitted. 

Report  Approved  by  _ 

M.  Duane  Horton 

Principal  Investigator 


R.  L.  Coates 
Co- Investigator 


;:Nt 


mcejsiw  Nr 

rmn  ibt 

Mff  SfCTIO*  □ 

HVMMnca  □ 

JtftTlflCUtM  . 

!*» 

OISTliiriM/ttllliiatlY  COOES  j 

OIST.  ,  AVAIL.  Ml.gr  tfECUll 

Qualified  Requestors  may  obtain  additional  copies  from  the 


/ 


Defense  Documentation  Center.  All  others  should  apply  to 


the  clearinghouse  for  Federal  Scientific  and  Technical  In¬ 


formation, 


TABLE  OF  CONTENTS 


Summary  . 

I.  Introduction  . 

II.  Radiation  Augmented  Burning  Rate  Studies  .  .  . 

1 .  Theory 

2 .  Experiment 

3.  Conclusions 

III.  Extinction  and  Other  Combustion  Characteristics 

1.  Oscillatory  Combustion 

2.  Lower  Deflagration  Permit 

IV.  Theory  . 

1.  Review  of  Previous  Models 

2.  Development  of  Present  Model 

V.  Comparison  Between  Theory  and  Experiment  .  .  . 

VI .  Conclusions  . 

VII.  References  . 

VIII.  Sample  Computer  Program  . 

IX.  Presentations  and  Publications  . 


ii 

1 

3 

7 

9 

20 

39 

40 
42 

52 


SUMMARY 

This  final  report  describes  the  progress  made  in  the  study  of 
propellant  extention  under  Grant  AF-AFOSR-0897-67  between  March  1967 
and  February  1969. 

At  the  inception  of  the  grant  period,  a  theory  had  been  developed 
to  describe  solid  propellant  flame  quenching  by  rapid  depressurization. 
Comparison  between  theory  and  experiment  seemed  promising,  and  it  was 
necessary  to  determine  the  magnitude  of  certain  combustion  parameters. 

An  experimental  technique  was  devised  and  used  to  determine  the  criti¬ 
cal  combustion  parameters.  The  use  of  these  parameters  then  revealed 
a  discrepancy  between  theoretical  predictions  and  experimental  results. 
In  fact,  the  discrepancy  was  so  great  as  to  show  the  theoretical  model 
to  be  invalid. 

Accordingly,  an  improved  theoretical  model  was  developed  and 
compared  to  experiment.  The  improved  theory  agreed  well  with  experi¬ 
ment  and  also  explained  the  causes  behind  some  previously  puzzling 
experimental  results. 


I .  INTRODUCTION 

This  report  describes  research  which  is  part  of  a  continuing 
program  intended  to  characterize  propellant  combustion  during  rapid 
depressurization.  The  research  was  initiated  under  Grant  AF-AFOSR-897-65 
and  continued  under  Grant  AF-AFOSR-897-67 . 

During  the  period  of  Grant  AF-AF0SR-897-cS,  a  simple  propellant 
extinction  theory  was  developed.  This  theory  was  based  on  a  combustion 
model  that  assumed: 

1 .  the  heat  conduction  is  one-dimensional 

2.  there  are  no  reactions  beneath  the  propellant  surface 

3.  the  thermal  properties  of  the  propellant  are  constant 

4.  the  propellant  is  homogenous  and  semi-infinite 

5.  the  heat  transfer  to  the  solid  propellant  surface  was  the  same 
as  the  steady  state  heat  transfer  at  the  same  pressure 

6.  the  surface  temperature  was  constant 

With  these  assumptions,  it  was  possible  to  numerically  solve  the  heat 
transfer  equations  and  describe  the  propellant  extinction  in  terms  of 
only  two  parameters.  A  dimensionless  heat  of  decomposition  and  char¬ 
acteristic  depressurization  time. 

Experimental  extinction  data  were  gathered  to  augment  the  data 
available  in  the  literature  and  then  these  data  were  compared  to  the  theory 
by  the  use  of  assumed  heats  of  vaporization.  The  comparison  was  favorable 
in  most  cases  but  depended  heavily  on  the  value  chosen  for  the  heat  of 
decomposition. 

A  simple  method  of  doubtful  reliability  was  used  to  experimentally 
determine  the  heats  of  decomposition  for  the  propellants  tested  in  this  program. 


Subsequent  use  of  these  values  to  make  theoretical  predictions  proved  discour¬ 
aging,  and  it  was  at  this  point  that  the  previous  grant  expired. 

It  was  then  necessary  to  find  a  better  method  for  determining  the  propel¬ 
lant  heat  of  decomposition.  This  was  accomplished  and  is  discussed  in  the 
next  section  of  the  paper.  Use  of  the  determined  heat  of  vaporization  showed 
that  the  theory  was  inadequate  and  hence  improved  extinguishment  theory  was 
derived.  This  theory,  based  on  the  Dennison  and  Baum  combustion  model  (1),  is 
described  in  section  IV.  As  is  shown  in  section  V,  the  improved  model  agrees 
well  with  experimental  observations. 
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II.  radiation  augmented  burning  rate  studies 


Since  the  exact  value  for  the  propellant  heat  of  decomposition  was  so 
critical  in  comparing  experiment  and  theory,  it  seemed  desirable  that  the 
parameter  be  determined  by  a  more  reliable  procedure.  Accordingly,  an  experi¬ 
ment.  was  devised  to  yield  data  from  which  the  burning  propellant  surface 
temperature  and  heat  of  decomposition  could  be  calculated. 

The  technique  involved  first  determining  the  burning  rate  of  propellant 
in  a  chamber  whose  wall  temperature  was  fi::ed  at  different  values  in  different 
tests.  Then,  the  burning  rate  of  the  propellant  was  determined  when  the 
propellant  was  the  shape  of  a  long,  thick-walled  cylinder  and  the  internal 
cylinderical  surface  burned.  It  was  then  reasoned  that  a  simple  strand  of 
propellant  burning  in  a  chamber  with  ccld  walls  experienced  a  net  radiative 
heat  loss  and  the  burning  rate  would  be  smaller  because  of  the  heat  loss. 
However,  as  the  chamber  wall  temperature  was  raised,  the  point  would  be  reached 
where  there  were  no  radiative  heat  losses  and  the  combustion  would  be  adiabatic. 
This  would  be  exactly  analogous  to  combustion  inside  a  long  hollow  cylinder 
where  the  radiative  losses  would  be  negligible.  The  propellant  surface 
temperature  could  be  calculated  from  the  chamber  wall  temperature  that  caused 
the  propellant  to  bum  at  the  same  rate  as  did  the  internally  burning  cylinder. 
For  the  PBAN-AP  propellant  tested,  Fig.  1  shows  the  surface  temperature  to  be 
560°C  in  excellent  agreement  with  the  experimental  data  of  others  (2,3). 

It  was  found  that  the  burning  rate  of  the  propellant  increased  linearly 
as  the  flux  from  the  chamber  wall  to  propellant  surface  increased.  A  computa¬ 
tional  scheme  was  developed  that  permitted  us  to  calculate  the  heat  of 
decomposition  for  the  burning  propellant.  The  scheme  is  outlined  as  follows: 
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FIGURE  1 


Tw4  CR4  x  |0-,z) 


-  Strand-burning  rates  as  a  function  of  wall  surface  temperature. 
The  dotted  line  corresponds  to  the  burning  rate  of  an  internal¬ 
burning  cylinder. 


1.  The  propellant  burning  rate  (r)  is  assumed  to  be  a  function  of  flame 
temperature  (T^)  and  radiant  flux  level  (f  )  as 

r  =  f(fr,Tf) 


2.  The  derivative  of  this  equation  is  found  to  be 


~ T~ 


r  *  (  STf  ) 


3.  A  heat  balance  is  made  at  the  interface;  the  gaseous  products  are  assumed 
to  have  a  characteristic  reaction  time  and  then  the  partial  derivatives 
are  written  in  terms  of  the  unknown  heat  of  decomposition  (y)  and 
several  parameters  whose  values  are  known. 

dTf 

4.  is  evaluated  from  an  overall  heat  balance, 
r 

5.  The  above  results  are  combined,  small  order  terms  neglected  and  the 
resultant  equation  is 


2PSICSCTS 


l  l  Cs(Ts  -  T0)  ♦  l] 

'  -  Tq)  TXT  L  1  +  Cg(Tfo  -  Ts)  J 


when  p  is  the  propellant  density,  C  is  its  heat  capacity,  C  the 
s  s  g 

mean  heat  capacity  of  the  gases  and  7  ,  TQ,  and  T^q  are  respectively 
the  surface  temperature,  the  conditioning  temperature,  and  the  adia¬ 
batic  flame  temperature. 

6.  The  value  for  the  heat  of  decomposition,  L,  is  found  by  rearranging 
the  equation  and  using  known  values  for  the  parameters  except  dr/df^, 
which  is  the  slope  of  the  heat  flux-burning  rate  curve  and  is  evalu¬ 
ated  from  Fig.  1. 
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For  the  propellant  tested,  the  heat  of  decomposition  was  found  to  be 
-130  cal/g  in  excellent  agreement  with  previous  work  (2).  This  number  means 
that  there  is  an  exothermic  reaction  at  or  near  the  propellant  surface. 

The  fact  that  the  decomposition  process  was  exothermic  proved  to  be  very 
significant.  Previously,  the  decomposition  was  assumed  to  be  endothermic 
and  all  calculations  were  made  on  that  basis.  When  the  correct  exothermic 
value  was  used  in  the  theory  and  the  theoretical  predictions  examined,  they 
were  found  to  be  totally  unrealistic.  In  fact,  it  was  found  that  a  decrease 
in  pressure  was  predicted  to  cause  an  increase  in  burning  rate.  A  relation 
total]''  at  variance  with  the  known  facts. 


III.  EXTINCTION  AND  OTHER  COMBUSTION  CHARACTERISTICS: 


A  readily  extinguishable  propellant  is  one  whose  burning  rate  is  easily 

made  to  approach  zero  by  a  given  pressure  transient.  Mathematically,  the 
dr/r 

value  of  is  large  for  such  a  propellant.  Similarly,  a  propellant  very 
susceptible  to  oscillatory  combustion  has  a  large  value  for  the  real  part  of 

dy/  p 

the  response  function  The  similarity  of  the  two  groups  suggests  that  a 

very  extinguishable  propellant  may  also  be  one  very  prone  to  oscillatory  com¬ 
bustion  . 

This  potential  relationship  was  investigated  by  determining  the  extinc¬ 
tion  characteristics  of  several  propellants  whose  oscillatory  combustion  was 
studied  elsewhere  (4,5).  Specifically,  a  mesa-type  base  (X-14)  a  conventional 
double  base  (JPN)  a  polyurethane  (JPL  534),  and  a  polysulphide  (T-35),  were 
tested.  Each  of  these  propellants  has  been  used  in  a  developmental  motor  and 
each  has  exhibited  severe  oscillatory  combination. 

The  experimental  results  are  shown  in  Table  1  and  clearly  indicate  that 
there  is  no  simple  relation  between  extinguishability  and  instability.  In 
fact,  the  two  double-base  propellants  are  about  equally  unstable  but  differ 
in  ease  of  extinguishability  by  two  orders  of  magnitude. 


Propellant 


( 


TABLE  1 
dr/r  V 

3p7p  ) 

/ 


Re 


Critical  d^p- 


JPN 

2.6 

-  7 

X-14 

2.4 

-  700 

JPL  534 

1.2 

-  20 

T-35 

1.4 

-  80 
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It  was  pointed  out  previously  (9)  that  there  is  a  correlation  between 
extinguishability  and  the  burning  rate  near  the  deflagration  limit.  However, 
the  burning  rate  near  the  deflagration  limit  is  a  parameter  rarely  measured 
and  reported.  Therefore,  during  this  period  little  additional  data  have  been 
found.  Such  data  as  were  found  supported  the  correlation  but  it  is  concluded 
that  at  this  time  the  use  of  such  a  seldom-measured  parameter  has  limited 
utility. 
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IV.  THEORY 


1.  Review  of  Previous  Models 

Mathematical  models  of  the  process  of  rapid  depressurization  extinguishment 
of  a  burning  solid  propellant  have  been  proposed  by  at  least  a  dozen  different 
authors.  According  to  a  recent  paper  by  Merkle,  et.  al.  (6),  the  universal 
approach  has  been  to  assume  that  the  extinguishment  process  occurs  because  of 
a  deficiency  in  the  heat  flow  from  the  flame  to  the  burning  surface  and  into 
the  unbumed  solid,  this  heat  flow  being  a  function  of  the  pressure.  This 
statement  applies  to  the  present  approach. 

To  describe  the  heat  flow,  the  solid  propellant  is  considered  to  be 
represented  as  a  semi-infinite  slab  in  order  that  only  one-dimensional  heat 
transfer  need  be  considered.  Figure  (2)  shows  the  geometrical  model  of  the 
process  along  with  the  general  shape  of  the  temperature  profile,  which  extends 
from  the  completely  reacted  products  of  combustion  at  the  adiabatic  flame 
temperature,  T^,  to  the  initial  conditioning  temperature  of  the  propellant, 

Tq,  For  further  simplification,  the  burning  surface  is  taken  as  the  reference 
point  for  the  position  coordinate  x,  the  propellant  assumed  to  move  to  the 
left  at  its  burning  rate  in  order  to  maintain  the  burning  surface  at  a  fixed 
position. 

An  energy  balance  taken  on  an  element  of  the  solid  of  thickness  dx  results 
in  the  following  partial  differential  equation: 


3T 

3t 


r 


3T 

dX 


CD 


where  a  is  the  effective  thermal  diffusivity  of  the  solid  and  r  is  the  burn¬ 
ing  rate.  This  equation  is  easily  integrated  fo:.*  the  steady-state  (3T/9t  =  0) 
to  yield  the  following  equation  for  the  temperature  profile  in  the  solid: 


(2) 


T  =  Tq  +  (Tg  -  T  )  exp  C -r  x  /  a) 

where  Tg  is  the  temperature  of  the  burning  surface  and  r  is  the  steady-state 
burning  rate.  The  heat  flux  that  must  be  supplied  to  the  solid  to  maintain 
this  profile,  in  addition  to  that  necessary  to  gasify  or  decompose  the  solid, 
is  given  by 

fs  *  "  V  CTs  -  V  (3) 

and  the  thermal  energy  stored  in  the  solid  above  its  initial  conditioning 
state  is 


qs  =  (k  /  r)  (Ts  -  TQ)  (4) 

where  k  is  the  thermal  conductivity. 

The  original  models  of  Von  Elbe  (7)  and  Paul  (8)  were  the  first  to  consider 
the  non-steady-state  heat  flow  problem  with  the  objective  of  predicting 
extinguishment  conditions.  These  models  are  very  similar  and  do  provide 
insight  into  the  physical  processes  that  might  lead  to  extinguishment  even 
though  they  are  based  on  rather  gross  simplifying  assumptions.  The  Paul  model 
is  outlined  here  to  provide  a  historical  reference  for  comparison  with  the 
model  developed  for  this  study. 

Paul  assumed  that  as  a  first  approximation  the  surface  temperature 
remained  constant.  He  further  took  the  steady-state  burning  rate  to  be  related 
to  the  pressure  by 

r  =  a  pn  (5) 

and  assumed  that  the  flux  into  the  solid  during  either  steady-state  or  transient 
periods  could  be  expressed  as 
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*,  -  *  VfTs  '  V 

For  a  transient  case,  t  would  be  interpreted  as  the  steady-state  rate  cor¬ 
responding  to  the  instantaneous  value  of  the  pressure.  During  a  transient 
period,  the  energy  stored  in  the  solid  must  change  as  the  burning  rate  changes, 
Both  Von  Elbe  and  Paul  assumed  as  a  first  approximation  that  this  stored 
energy  adjusts  instantaneously  with  changes  m  burning  rate,  its  magnitu  e 
always  being  given  by  equation  4.  Thus,  during  a  transient  period,  the 

required  flux  is 

f$  =  P  csr(Ts  -  Tq)  +  dqs/dt  (7) 

Combining  (4)  and  (5)  and  carrying  out  the  differentiation  yields 

.  /At  _  —  fT  -  T  )  (8) 

dqs/dt  -  -  r-  V  dt 


and  combining  (6) ,  (7) ,  and  (8)  yields 


r  -  r 


Equation  9  predicts  that 


[1  .  CWr*5  (9) 

the  instantaneous  rate  will  drop  to  zero  whenever 


(  =  -  r/na 

\  dt  / Ext. 


Since  it  was  observed  that  this  equation  only  qualitatively  correlated 
experimental  extinguishment  data,  a  correction  factor,  X,  was  added  to  bring 
it  into  quantitative  agreement,  or 


/«!L2 1 

V  dt  /Ext. 


=  -  Ar/nct 


a 


Horton  (9)  recognized  that  a  serious  shortcoming  of  the  Paul  and  Von  Elbe 
theories  was  the  assumption  that  the  heat  stored  in  the  solid  adjusted 
instantaneously  to  changes  in  burning  rate.  In  reality  when  a  change  in 
pressure  and  corresponding  heat  flux  to  the  solid  occurs  there  is  a  time  lag 
before  a  new  steady-state  profile  is  established.  To  treat  the  transient 
conduction  process  more  exactly,  Horton  carried  out  a  numerical  integration 
of  the  non-steady-state  differential  energy  balance  equation  (Equation  1) . 

He  further  improved  the  theory  by  including  the  energy  required  for  gasifica¬ 
tion  of  the  solid  in  the  solid  heat  flux  equation.  His  equation,  corresponding 
to  Equation  6,  for  the  flux  supplied  to  the  solid  is 

f s  =  r  p  [L  *•  Cs(Ts  -  Tq)]  -r  pL  (12) 

where  L  is  the  effective  heat  of  gasification  of  the  solid.  The  flux  required 
by  the  solid  in  Horton’s  model,  corresponding  to  Equation  7,  is  obtained  from 
the  numerical  integration  of  the  differential  energy  balance.  As  in  the  Paul 
model,  the  criterion  for  extinguishment  is  that  the  rate  becomes  zero  as  a 
consequence  of  the  supplied  fg  being  entirely  absorbed  by  the  solid  in  chang¬ 
ing  qg,  the  amount  of  heat  stored. 

Horton's  model  predicts  for  an  exponential  pressure  decay  and  for  a 

given  and  Tg  that 

(d  InjA  .  .  wP/cx)  (13) 

\  dt  /Ext.  H^ 

where  AH  is  a  function  of  n,  Lg,  and  Pa/PQ,  the  latter  being  the  ratio  of  the 
final  pressure  following  depressurization  to  the  initial  pressure.  The  function 
v  is  described  in  terms  of  graphs  developed  from  the  results  of  the  numerical 


calculations. 


Comparing  (13)  with  (11)  shows  that  X^  corresponds  to  the  correction 
factor  X  employed  with  the  Paul  theory.  Thus,  Horton's  theory  provides  a 
method  of  evaluating  X  in  ter, .vs  of  more  fundamental  parameters. 

During  the  course  of  the  present  contract,  it  was  discovered  that  Horton's 
model  yielded  erroneous  results  if  the  heat  of  gasification,  L,  were  negative. 
This  was  considered  to  be  a  serious  shortcoming  since,  as  discussed  in 
Section  II,  there  is  evidence  suggesting  this  parameter  can  indeed  be  negative 
for  some  propellants. 

More  recent  models  have  been  proposed  by  Wooldridge  (10)  and  by 
Summerfield  (11).  The  latter  model  is  based  upon  the  Granular  Diffusion  Flame 
Theory  of  Summerfield,  which  does  not  provide  the  proper  dependence  of  steady- 
state  burning  rate  on  pressure  for  many  propellants.  Wooldridge's  approach  was 
to  modify  a  combustion  theory  originally  proposed  by  Denison  and  Baum  (1). 

This  same  approach  has  been  followed  in  the  present  study.  The  unique  feature 
of  the  present  approach  is  that  the  Denison-3*v<»:<!  theory  is  used  only  to  predict 
deviations  from  the  steady-state  conditions;  the  conventional  steady-state 
strand -burner  data  is  employed  as  fundamental  input  data. 

2.  Development  of  Present  Model 

The  Denison-Baum  model  was  originally  derived  as  a  simplified  transient 
combustion  model  for  application  to  combustion  instability  problems.  Rather 
than  employing  a  constant  burning  surface  temperature,  as  Horton  did,  this 
model  uses  an  Arrhenius-type  relationship  between  burning  rate  and  temperature. 

r  =  A  exp  (-E  /RT^)  (14) 

It  is  further  assumed  that  the  rate  is  coupled  to  the  flame  temperature, 
Tp  and  pressure  by  the  following  equation 
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exp  (Ef/RTf) 


(15) 


T  N+l 
f 


In  the  present  analysis  these  equations  have  been  replaced  with 

r  =  rp  exp  [(-Eg/R)  (1/Ts-1/Ts,  )]  (16) 

and 

r  =  rp(T£/Tf)n+1  exp  [ (-Ef/R) (1/Tf  -  1/T£,  )]  (17) 

where  rp,  Tg,p,  and  T£,p  are  the  steady-state  burning-rate,  surface  temperature, 
and  flame  temperature  at  the  instantaneous  pressure  p,  and  n  is  the  empirical 
value  of  dlnrp/dlnp. 

The  most  important  feature  of  the  Denison-Baum  theory  is  the  equation  for 
the  rate  of  heat  conduction  into  the  solid  beneath  the  burning  surface.  This 
equation  may  be  written  as 

k  f|x  = 0  *  =  -r°stcs(Ts  -  v  *  Vrf  -  Vi  tl8) 

with  the  origin  of  the  x  coordinate  fixed  on  the  burning  surface,  the  positive 
direction  extending  into  the  solid.  This  equation  is  obtained  by  combining 
overall  energy  balances  for  both  steady  and  non-steady  conditions,  eliminating 
energy  terms  associated  with  the  surface  gasification  reaction.  The  virtue  of 
this  equation  is  that  data  for  surface  and  flame  temperatures  are  more  readily 
available  than  data  for  effective  heats  of  gasification.  Calculation  of  the 
numerical  value  of  the  instantaneous  heat-flux  supplied  to  the  solid  is  also 
straight-forward.  Given  instantaneous  values  of  Tg  and  P,  r  is  computed  from 
(16)  and  T£  is  computed  from  (17) .  These  values  then  given  the  instantaneous 
value  of  the  flux  to  the  solid  when  substituted  into  (18) . 
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The  most  difficult  problem  in  employing  either  the  Horton  or  the  Denison- 
Baum  models  in  predicting  extinguishment  conditions  is  the  calculation  of  the 
rate  of  transient  heat  conduction  into  the  solid  propellant.  The  usual 
numerical  technique  for  solving  the  differential  energy  balance  which  describes 
the  transient  heat  transfer  process,  Equation  (1),  is  to  separate  the  solid 
into  finite  difference  elements  such  that  the  energy  balance  on  each  element 
can  be  written  as  an  algebraic  equation  rather  than  a  differential  equation. 

The  algebraic  equations,  one  for  each  of  the  finite  elements,  are  then  solved 
simultaneously  to  obtain  both  the  temperature  profile  at  any  given  instant 
during  the  transient  period  2nd  the  increment  of  change  at  any  position  for  a 
given  increment  in  time. 

Since  sophisticated  numerical  techniques  have  become  available  for  the 
simultaneous  solution  of  a  set  of  ordinary  differential  equations,  two  recent 
authors  (12,  13)  have  pointed  out  the  advantage  of  using  these  techniques 
to  solve  transient  heat  conduction  problems.  The  differential  energy  balance 
for  both  finite  elements  is  transformed  into  an  ordinary  differential  equation 
applying  finite  difference  methods  to  the  space  derivatives  only.  Available 
differential  equation-solving  programs  are  then  employed  to  carry  out  simul¬ 
taneous  integration  of  these  equations  to  compute  the  variation  of  the  tempera¬ 
ture  profile  with  time.  This  approach  has  been  adopted  in  this  study. 

Figure  2  illustrates  the  method  of  separating  the  solid  into  finite 
elements.  The  ordinary  differential  equation  which  approximates  the  energy 
balance  for  each  element  is  written  as 


dTi 

~dt 


«t- 


Ti-i  -  ni +  Ti.i,  .  ,Ti.i  -  h-i, 

— ]  *  rt — m — 'i 


(AX) 


(19) 
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The  wide  variation  of  r  during  the  extinguishment  transient  presents 
special  problems  in  the  numerical  solution  of  the  transient  conduction  equa¬ 
tion.  This  problem  is  illustrated  by  considering  the  change  in  the  depth  of 
penetration  of  the  temperature  profile  during  the  transition  from  one  steady- 
state  condition  to  another,  the  first  condition  being  that  prior  to  depressuri¬ 
zation  and  the  other  the  condition  after  a  depressurization  has  occurred,  the 
propellant  assumed  to  continue  burning.  Tie  steady-state  temperature-profile 
is  described  by 

f  - 

7f — -y~Y~  =  exP  (-rx/°0  (20) 

s  03 

Suppose  that  during  the  transient  the  rate  is  reduced  by  a  factor  of  10. 

Equation  (21)  shows  that  the  depth  at  which  a  given  temperature  is  established 
in  the  solid  is  then  increased  by  a  factor  of  10.  Suppose  further  that  20 
equal  elements  of  width  AX  are  selected  at  the  start,  the  final  element  chosen 
such  that  its  temperature  rise  was  only  5%  of  the  rise  at  the  surface.  The 
temperature  rise  of  this  element  when  the  new  steady-state  is  established  would 
be  74%  of  the  surface  rise;  200  elements  of  width  AX  would  now  be  required  to 
reach  a  depth  where  the  rise  was  only  5%.  Since  a  separate  equation  must  be 
solved  for  each  element,  the  problem  can  become  unwieldy  if  equal-sized  elements 
are  employed  and  the  burning  rate  change  is  large. 

Marxman  and  Wooldridge  (10)  suggests  that  the  problem  of  excessive  finite 
elements  can  be  circumvented  by  employing  the  transformation 

1  -  Y  =  exp  (-rx/a)  (21) 

to  define  a  new  distance  variable  Y.  This  variable  ranges  from  0  to  1  as  X 
ranges  from  0  to  »,  and  elements  of  equal  width  AY  correspond  to  AX  elements 
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which  become  successively  wider  as  X  increases,  the  last  element  having  the 
width  This  transformation  has  been  investigated  during  this  study  and  found 
to  lead  to  erroneous  results  if  the  rate  varies  widely.  A  different  trans¬ 
formation  has  been  found  to  better  accomplish  the  desired  reduction  in  finite 
elements  and  also  retain  the  accuracy  of  the  numerical  integration  method.  In 
this  transformation  the  new  distance  variable,  y,  is  defined  by 


P-  =  (1  -  y) 

dx  a 


(22) 


The  corresponding  finite  element  spacing  is 

Ax  =  (a/r)  Ay/(1  -  y)  (23) 

Letting  y  range  between  zero  and  one  and  taking  Ay  as  a  fixed  quantity,  the 
finite  difference  spacing.  Ax,  increases  regularly  with  both  increasing  depth 
beneath  the  surface  and  decreasing  burning  rates.  The  resulting  finite 
difference  equation  in  non-dimensional  form  that  was  used  in  place  of  Equation  (1) 
is 


aiR2(ei-rV9i.i)-bici(ei.rei-i 


dR 

)  (1/R)  3f 


(24) 


with 


R  +  r/r 
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(25) 


ai  =  (1  -  y)0  -  yi-Ay/(Ay)2 

=  1  Ay/(1  -  Y:) 

1  .  z  1 

3=0 

C.  =  (1  -  yi)(l  -  y.  -  Ay)/(2  -  2y.  +  Ay) 

The  last  term  on  the  right-hand  side  of  Equation  (25)  accounts  for  adjustment 
of  the  temperature  as  the  mid-point  of  the  finite  element  is  shifted  along  the 
x  axis  due  to  changes  in  r.  This  term  does  not  appear  in  the  equation  for  the 
surface  temperature  (which  becomes  Tp  since  this  point  is  stationary. 

The  surface  boundary  condition.  Equation  (18),  was  expressed  in  finite- 
difference  form  by  defining  a  station  at  a  position  y  above  the  surface  with 
a  fictitious  temperature  Tq: 

T  =  T  +  2  Ax  (kS  (26) 

0  2  0  3x  x=o+ 

where  is  the  temperature  at  the  first  station  beneath  the  surface. 

Initial  calculations  were  made  solving  Equations  16,  17,  18,  24,  25, 
and  26  simultaneously  along  with  the  following  equation  describing  an  imposed 
depressurization  rate. 


f-=  -  (P  -  P  }  (0.693/-  )  C27) 

s 

A  standard  computer  program  for  integrating  simultaneous  ordinary  differential 
equations  based  on  Adams  method  (15)  was  employed  The  Fortran  Coding  for 
carrying  out  these  calculations  using  the  BYU  Librascope  L-3055  computer  is 
included  in  the  Appendix. 


V.  COMPARISON  BETWEEN  THEORY  AND  EXPERIMENT 

A  series  of  parametric  calculations  have  been  made  both  to  determine  a 
consistent  set  of  parameters  for  correlating  observed  extinguishment  conditions 
and  to  determine  the  sensitivity  of  the  model  to  variations  in  the  individual 
parameters . 

A  rough  correlation  of  existing  depressurization  extinguishment  data  is 
shown  in  Figure  3,  where  the  time  corresponding  to  reducing  the  pressure  to 
one -half  the  initial  value  is  plotted  versus  the  initial  burning  rate.  The 
data  po  nts  represent  limiting  conditions  for  extinguishment  to  occur.  If  t^ 
is  greater  than  that  indicated  by  the  data  point  and  the  initial  rate  is  the 
same,  the  propellant  will  continue  to  burn  rather  than  extinguish. 

The  line  correlating  this  data  can  be  represented  by  the  equation 

tj^  =  0.5 /F2  (28) 

where  t,__  is  in  milliseconds  and  F  is  in  inches  per  second.  The  rate  of 
depressurization  is  assumed  exponential  with  P  /P  assumed  negligible  compared 
to  unity. 

The  non-dimensional  half-time  corresponding  to  (28)  is 

=  0.5(10"3)/a  (29) 

Thus,  if  a  is  again  assumed  to  be  2.5(10~4)  in.2/sec.,  the  value  of 
corresponding  to  extinguishment  is  2. 

Figure  4  shows  predicted  burning  rates  during  depressurization  transients 
for  Tj  =1.5  using  the  following  set  of  parameters: 

Es  =  30  Kcal/gm.  mole 

Eg  =  30  Kcal/gm.  mole 
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Initial  Burning  Rate  -  in/sec 

FIGURE  3  -  Results  of  attempt  to  correlate  marginal  p 
extinguishment  with  initial  burning  rate 


T  =  298°K 
o 

T,  =  10  T 
t  o 

T  =  3  T 
s  o 

n  =  0.4 

Pa/Po  =  0.1,  0.01 

Cg/c  =  1.0 

The  burning  rate  is  shown  to  establish  a  new  steady-state  value  following 
depressurization  to  P  /P  =0.1;  however,  it  is  shown  to  continue  to  decrease 

ci  O 

following  depressurization  to  Pa/PQ  =  0.01,  the  ratio  r/F  dropping  far  below 
the  value  corresponding  to  steady-state  at  the  pressure  following  depressuriza¬ 
tion.  This  latter  behavior  is  considered  to  represent  the  extinguishment 
process,  the  reasoning  being  that  combustion  will  not  be  sustained  below  some 
absolute  minimum  rate  (or  absolute  minimum  Ts).  For  example,  if  the  initial 
rate  were  0.2  inches  per  second,  the  predicted  rate  for  the  case  of  P&/ PQ  = 

0.01  at  t  =  24  is  0.0014  inches  per  second.  Since  rates  this  low  are  not 
observed  experimentally  during  steady-state  conditions,  it  seems  improbable 
that  rates  this  low  would  exist  during  a  transient  condition.  Thus  it  is 
predicted  that  if  =  1.5  extinguishment  would  become  marginal  at  some  value 
of  Pa/PQ  lying  between  0.01  and  0.1  for  this  particular  set  of  propellant 
parameters,  and  since  =  2.0  is  required  to  correlate  the  experimental  data, 
as  discussed  above,  this  particular  combination  of  parameters  yields  reasonably 
good  agreement  between  theory  and  experiment. 

A  series  of  computations  has  been  made  to  determine  the  effect  of  varying 
the  gas-phase  and  solid-phase  global  activation  energies.  The  results  of  these 
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calculations  are  shown  in  Figures  S  and  6.  Increasing  either  Eg  or  is 

shown  to  result  in  greater  susceptibility  to  extinguishment.  This  effect  is 

most  strikingly  shown  in  Figure  5  where  with  E  =  15  a  new  steady-state  burning 

8 

rate  is  established  following  depressurization  while  with  E^  =  30  the  rate 
continues  to  decrease. 

Figures  7,  8,  9,  and  10  present  the  results  of  a  series  of  calculations 
where  P  /P  and  x,  were  systematically  varied.  These  calculations  show  a 

cl  O 

rather  strong  dependency  of  the  value  of  x^  required  for  extinguishment  on 
P  /P  .  Assuming  that  extinguishment  corresponds  to  the  case  where  r/r  is 
predicted  to  drop  below  0.05,  the  following  table  summarizes  the  predicted 
conditions  for  extinguishment  derived  from  this  set  of  calculations. 


Pa/PQ  (t^)  Extinguish 


0.01 

10-20 

0.05 

1-3 

0.10 

0.1-0. 5 

0.50 

<  0.01 

These  results  suggest  that  the  X!  required  for  extinguishment  varies  roughly 
with  the  inverse  square  of  Pa/PQ .  This  information  may  provide  some  explanation 
for  the  scatter  in  the  correlation  of  experimental  data  shown  by  Figure  3. 

The  model  for  the  solid  propellant  combustion  process  has  also  been 
combined  with  equations  describing  the  transient  ballistics  of  a  solid  rocket 
motor.  Equations  expressing  conservation  of  both  mass  and  energy  inside  the 
motor  cavity  have  been  derived  previously  (14).  In  the  present  study  these 
equations  have  been  used  in  the  following  form: 
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FIGURE  8 


Predicted  Effect  of  Varying  tj.  for  Pa/P 


•Wf 


0  4  8  12 

Dimensionless  Time  -  t  r2/o C 


FIGURE  9  -  Predicted  Effecl  of  Varying  for  Pa/PQ  =  0.1 
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(30) 

(31) 

with 

P  =  p/p 


dP/dxc  =  YR9f/0c)  (l-(9c/6f)  (Rjj/R) 

_  YOf 

d0c/dxc  =  C0c/pH0c/ec)(R)[i-Y+(ri  -DfR/V1 


A  ■  <Avent  *  V/At 

Prediction  of  the  instantaneous  pressure  versus  time  requires  the 
simultaneous  numerical  integration  of  these  equations  along  with  those 
discussed  earlier  for  the  solid  and  the  use  of  the  steady-state  burning 
relation.  Provision  has  been  made  in  the  program  to  use  the  Summerfield 
Equation,  the  Vielle  equation,  or  a  table  relating  the  burning  rate  to 
pressure.  In  all  three  cases,  extinguishment  is  taken  to  occur  if  the  rate 
becomes  less  than  the  burning  rate  at  the  deflagration  limit  of  the  propellant. 

Tables  II  and  III  show  results  of  calculations  made  with  the  combined 
combustion-ballistics  model. 
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Table  II 


Predicted  effect  of  various  parameters  on  extinguishment  of  a  propellant  that 
follows  the  Summerfield  burning  rate  equation 


Parameter  Varied 


Conditions  for  Marginal  Extinguishment 
Vent -Area  Rates  Initial  Depressurization  Rate 

VAi  -Cp/p)  o 


Reference* 

2. 8-2.9 

145-152 

To  =  100°C 

3. 6-3. 7 

204-211 

E  =  E  =25 
s  g 

2. 7-2. 8 

137-145 

E  =  28 
s 

2. 7-2. 8 

137-145 

E  =  28 
g 

2. 7-2. 8 

137-145 

E  =  35 
s 

2. 7-2. 8 

137-145 

E  =  35 
g 

2. 5-2. 6 

125-130 

T  =  382 
s 

3. 7-3. 8 

211-218 

Tf  =  2185°K, 

L*  =  364,  C*  =  3910 

2. 6-2. 7 

101-107 

Tf  =  2185°K, 

C* 

'  =  3910,  Kn  =  274 

2. 5-2. 6 

108-115 

Kn  =  275,  L* 

= 

364,  Po  =  667 

3. 3-3.4 

159-166 

Kn  =  325,  L* 

a 

432,  Po  =  886 

3. 9-4.0 

171-176 

Kn  =  188,  L* 

a 

250,  Po  =  340 

2. 2-2. 3 

121-131 

Kn  =  375,  L* 

= 

497,  Po  =  1124 

4. 5-4. 6 

179-184 

Kn  =  325,  Po 

= 

886 

3. 6-3. 8 

203-218 

Kn  =  375,  Po 

a 

1124 

4. 0-4.1 

232-239 

Kn  =  188,  Po 

= 

340 

2. 4-2. 5 

106-114 

Kn  =  275,  Po 

= 

667 

3. 2-3. 3 

174-182 

L*  =  1000 

4. 0-4.1 

80-83 

L*  =  105 

2. 2-2. 3 

254-274 

C  /C  =  1.0 
s  g 

3. 2-3. 3 

174-182 

Kn  =  274,  C* 

a 

3910 

2. 5-2. 6 

108-115 

*The  reference 

To=25°C 

Kn=238 

L*=316  in 

Tf=2910°K 

C  /C  =.75 
s  g 

Pq=520  psia 

r  =.25  in/sec 
o 


propellant  and  motor 

T  =600°C 
s 

Ps=.062  lbs/in3 

Y=1 . 22 

C*=4S00  ft/sec 

E  =21  K  cal 
s 

E  =21  K  cal 
g 


had  the  following  properties: 

ambient  pressure=15  psia 

propellant  area=8.3  in2 

nozzle  open  time=.001  sec 

burning  rate  at  1000  psi=.33  in/sec 

exponent  at  1000  psi=.4 

burning  rate  at  deflagration  limit= 

.009  in/sec 
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Table  III 


Effect  of  various  parameters  on  extinguishment  of  a  propellant  that  follows 
the  Vielle  burning  rate  equation 


Parameters  different  than  reference 


Extinguishment 

VA1 


Occurs  between 
“(P/P)n 


Reference* 
n  =  .7 
n  =  .3 

burning  rate  at  deflagration  limit  =  .003 


Kn  =  415> 

L*  =  582 

ro  =  .144 

Kn  =  166, 

L*  =  221 

ro  =  .360 

K  =  1190 
n 

,  L*  =  1580 

ro  =  .0495 

K  =  131, 
n  ’ 

L*  =  174 

ro  =  .480 

3. 2-3. 3 

1.7- 1. 8 
9.0-9.25 

3. 2- 3. 3 

2. 7- 2. 8 

3. 7- 3. 8 

2. 3- 2.4 
4. 2-4. 3 


177-185 

59-68 

554-569 

177-185 

78-82 

292-303 

23-24 

425-499 


*The  reference  propellant  and  motor  had  the  following  properties: 


To  =  25°C 

K  =  238 
n 

L*  =  316  in 

Tf  =  2910°K 

C  /C  =  .75 
s  g 

Pq  =  520  psia 

r  =  .25  in/sec 
o 


T  =  600°C 
s 

Ps  =  .062  lbs/in3 
Y  =  1.22 

C*  =  4500  ft/sec 

E  =  21  K  cal 
s 

E  =  21  K  cal 
g 


ambient  pressure  =  15  psia 
propellant  area  =  8.3  in2 
nozzle  open  time  =  .001  sec 
rate  at  520  psi  -  .25  in/sec 
n  =  .5 

r  at  deflagration  limit  =  .009 
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One  of  the  parameters  noted  to  have  a  significant  influence  on  marginal 
extinguishment  conditions  was  the  initial  L*  of  the  motor.  Additional  calcula¬ 
tions  in  which  this  parameter  was  varied  by  a  factor  of  4  again  showed  the 
importance  of  the  overall  pressure  drop.  In  these  calculations  the  initial 
burning  area/throat  area  ratio  and  initial  depressurization  rate  were  held 
constant  by  adjusting  the  venting  area.  The  results  of  these  calculations 
are  shown  in  Figure  10.  Because  a  much  larger  vent  area/ throat  area  ratio  is 
necessary  in  the  larger  L*  motor  to  achieve  the  same  initial  depressurization 
rate  as  the  smaller  L*  motor,  the  final  pressures  differ  in  the  two  cases  as 
do  the  final  portion  of  the  blow-down  curves.  The  larger  area  ratio  is  seen 
to  be  sufficient  to  cause  the  pressure  in  the  larger  motor  to  drop  to  near 
ambient,  the  propellant  burning  rate  dropping  to  a  level  consistent  with 
extinguishment.  The  pressure  in  the  smaller  L*  motor  levels  off  at  a  new 
steady-state  well  above  ambient  with  no  indication  that  extinguishment 
would  occur.  If  the  model  predicts  a  response  that  is  truly  indicative  of 
the  actual  response  in  a  motor,  as  suggested  by  the  comparisons  made  above, 
this  set  of  calculations  clearly  shows  the  difficulty  that  might  be  encountered 
in  correlating  data  for  extinguishment  conditions  described  only  in  terms  of 
the  initial  depressurization  rate  or  t^.  Critical  initial  rates  obtained  with 
end-burning  motors  which  usually  have  large  L*’s  should  not  be  expected  to  be 
the  same  as  critical  initial  rates  obtained  with  motors  containing  internal 
burning  cylinders,  which  normally  have  low  L*'s. 

Some  data  from  References  10  and  12  are  shown  in  Figure  11  and  tend  to 
confirm  the  predictions  that  have  been  made.  That  is,  the  initial  depressurization 
raze  required  to  extinguish  a  propellant  in  a  given  motor  may  be  significantly 
different  from  that  in  a  different  motor--thus  making  design  on  the  basis  of 
dP/dt.  data  very  difficult. 
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P/P  sec- 


PU  269  UTX  10645 

End  Burning  Motor  □  1 5  Swing  Nozzlo  Motor  0 

Tubular  Motor  ■  5*  Swing  Nozzlo  Motor  O 

Doublo  Brain  Window  Motor  0 
Window  Motor  0 

End  Burning  Window  Motor  0 
Slab  Motor  0 


r0  in/sec 


Data  taken  from  Reference*  10  and  12  showing  that  propellant  In  different 
motors  requires  differenct  depressurization  rates  to  be  extinguished. 


Figure  11 


Figure  12  shows  that  experimental  data  for  the  effect  of  motor  size  can 
be  correlated  and  also  that  the  predictions  of  the  improved  extinguishment 
model  are  in  good  agreement  with  the  experimental  data.  However,  Figure  12 
probably  does  not  represent  the  best  form  of  correlation  that  can  be  made. 

As  is  discussed  in  Reference  16,  it  appears  that  a  simpler  and  better  correla¬ 
tion  relates  the  final  L*  of  the  motor  to  the  final  burning  rate.  That  is, 
the  steady-state  burning  rate  that  would  result  from  the  increased  nozzle  area 
if  stable  combustion  resulted.  Figure  13  shows  results  of  the  parametric 
study  presented  in  this  type  of  correlation. 


-35- 


Vielle  Equation  O 
Summerfield  Equation  ■ 


4  6  8  .10 


4  6  8  1.0 


rf  (in/sec) 

Data  taken  front  Reference  (10)  and  (12)  showing  a  relation  between  initial  pressure 
(P  ) ,  depressurization  rate  (Po) ,  burning  rate  (  ) ,  thermal  dif fusivity  (  ) , 

chamber  free  volume  (v) ,  and  exposed  propellant  surface.  (Ap) 

Figure  12 
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Theoretical  relationship  between  the  final  conditions  in  a  motor  that  marginally  fails 
to  be  extinguished,  for  two  types  of  steady-state  burning  behavior. 


Figure  13 


CONCLUSIONS 


1.  It  has  been  shown  that  the  effect  of  radiation  on  the  propellant  burning 
rate  can  be  used  to  evaluate  some  combustion  parameters  and  the  values  can 
then  be  used  to  test  the  consistancy  of  a  combustion  model. 

2.  Extinguishment  and  oscillatory  combustion  were  shown  not  to  be  related-- 
at  least  not  in  a  simple  manner. 

3.  An  improved  extinguishment  model  has  been  developed  and  shown  to  agree 
well  with  experimental  data. 
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